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ABSTRACT 

■ We review our recent results on the derivation of a B-S equation in QCD in a 

Wilson loop context. We work in a second order formalism, use the Feynman- 
Schwinger path integral representation for a quark in external field and obtain 
a similar expression for a quark-antiquark amplitude in which the gauge fields 
Q ■ appears only through the Wilson loop integral W. A qq B-S equation is obtained 

| starting from this expression under the assumption, already used in the derivation 

of heavy quark potential, that i In W can be expressed as the sum of a perturbative 
contribution and an area term. The intrinsically nonperturbative derivation method 
is first discussed on the simpler case of two spinless particles interacting through a 
J> ' scalar field. The qq semirelativistic potential is reobtained in the large quark mass 

limit. 
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Introduction 



Many attempts have been made to apply the Bethe-Salpeter equation to a study 
of the spectrum and the properties of the mesons, hoping to obtain an unified and 
consistent description of the quark-antiquark bound states envolving light quarks 
as well as heavy ones. In all these attempts, to our knowledge, the choice of the 
confinement part of the kernel was purely conjectural and only made in such a way 
that the successful heavy quarks potential could be recovered in the non relativistic 
limit. Beside being theoretically unsatisfactory, as well known, such a procedure 
encounters many conceptual and phenomenological difficulties 0. 

On the contrary, in the case of large quark masses it was possible to derive the 
semirelativistic qq potential ''essentially from first principles" reducing it to an eval- 
uation of the Wilson loop integral ™ 

W = ^(TrP T expig{f dafAn}) (1) 

(the loop T is supposed made by a quark world line (r\), an antiquark world line 
(T2) followed in the reverse direction and closed by two straight lines connecting the 
initial (2/1,2/2) and the final (xi,x 2 ) points of the two world lines; A^(x) denotes a 
colour matrix of the form A^x) = |A^A a ; Pr prescribes the ordering along the loop 
and Tr denotes the trace on the colour matrices; the expectation value stands for a 
functional integration on the gauge field alone (see below)). 

The basic tool for this derivation was the use of a path integral representation of 
the quark propagator in an external field and the basic assumption for the evaluation 



of W was that the logarithm of this quantity could be written as the sum of a 
perturbative contribution and an area law term 

ilnW = i(ln WOpert + <rS min . (2) 

In principle S m i n should be the minimum area enclosed by T; in practice one uses its 
straight line approximation, consisting in replacing it with the surface spanned by the 
straight lines connecting equal time points on the quark and the antiquark worldlines. 
Precisely one writes 

S mi n= dtr / 1 dA[l-(Av 1T + (l-A)v 2T ) 2 ]^, (3) 
Jti Jo 

(t stands for the ordinary time, Zi(t) and z 2 (i) for the quark and the antiquark 
positions at the time t, Vj = -Jjr denotes the ordinary velocity and the suffix T 

specifies its transverse component (5 hk — ^-^-)vj, r(£) denotes the relative position of 
the two particles zi(t) — z 2 (t)). Actually this equation can be shown to be exact up 
to the v 2 terms. 

As well known, due to asymptotic freedom, the first term in (2) is assumed to 
describe correctly the short range interaction and vanishes as the average distance f 
between the quark and the antiquark increases; the second one can be justified on 
the basis of the strong coupling expansion, is assumed to describe the long range 
interaction and vanishes as f decreases. Possibly (2) is too poor at intermediate 
distances and instanton contributions should be taken into account (attempts already 
exist in this direction, see e.g. H). Furthermore, presently the theory is not able to 
give a relation betweeen the "string tension" a and the renormalized strong coupling 

2 

constant a s = j-, which in practice have to be treated as independent variables. In 
spite of the above circumstance, significant results can be already produced by (2) 
and (3). 

The difficulty in obtaining a confining Bethe-Salpeter equation in QCD comes 
from the fact that usually this equation is derived by appropriate resummations of 
the perturbative series expansion and obviously such a procedure can provide only the 
perturbative part of the kernel. In this paper we want to review some of our recent 
results a to show that the technique used for the potential can be adapted even to the 
case of the Bethe-Salpeter equation and that a theoretically well founded kernel can 
be obtained in this way on the basis of (2) and (3). Since (3) is not Lorentz invariant, 
such equation shall be assumed to be valid in the center of mass frame. 

Let us begin to present our results in more precise terms. The QCD lagrangian 
can be written as 

L = Y, hd-fD^ - m f W f - -F, U F^ + L GF (4) 
/=i 4 

(where = — igA^ and Lqf is the gauge fixing term), and as usual the gauge 
invariant quark-antiquark Green function as 

1 — — c 

Gf(x 1 ,x 2 ,y 1 ,y 2 ) = -(0\T4) r 2 (x 2 )U(x 2 ,x 1 )ij 1 (x 1 )4) 1 (y 1 )U(y 1 ,y2)iJ 2 (y2)\0) = 



= ^Tr{U(x 2: Xl )SW( Xl , Vl ; A)U(y u y 2 )C- 1 S^(y 2: x 2 ; A)C)(5) 

Here c denotes the charge-conjugate fields, C is the charge-conjugation matrix, U the 
path-ordered gauge string 



U(b, a) = P ba exp (ig £ dx» A^x)j (6) 

(the integration path being along the straight line joining a to b), S\ and S 2 the quark 
propagators in the external gauge field A^ and obviously 

fV[A]M f (A)f[A]e iS W 
{J[ l) fV[A]M f (A)e iS W ' {l) 

S[A] being the pure gauge field action and Mf(A) the determinant resulting from the 
explicit integration on the fermionic fields (which however in practice we set equal to 
1 in the adopted approximation). 

The propagators Si and S 2 are supposed to be defined by the equation 

(i^D,-m)S(x,y;A) = 5\x-y) (8) 

and the appropriate boundary conditions. 

In the derivation of the potential one performs over (8) a Foldy-Wouthuysen 
transformation, replaces the Dirac type propagators S''- 1 -' and S^ by Pauli ones 
and takes advantage of ordinary path-integral representations for and 

K<®. In this way one obtains a two-particle Pauli propagator in which the gauge 
field appears only through the quantity W and its functional derivatives with respect 
to the world lines of the quark q and the antiquark q. Finally, after using (2) and (3) 
and expanding in the velocities one can verify that the obtained propagator satisfies 
a two-particle Schrodinger equation. 

For a fully relativistic treatment it is convenient first to pass to the second order 
formalism by setting 

S(x,y;A) = (i'fD v + m)A a (x,y;A), (9) 

with 

(D„D» + m 2 - l -g a^F^)A a (x, y; A) = -5\x - y) (10) 
(cr^ = |[7 M ,7^]). Then, taking into account gauge invariance, we can write 

Gf(xi,x 2 ;yi,y 2 ) = (r^d^ + mi)(ij 2 d X2U + m 2 )H A {xi, x 2 ; y u y 2 ) (11) 

with 

H 4 (xi,x 2 ;yi,y 2 ) = ~Tr(U{x 2 , xi)A^(x u Vl ; A)U(y u y 2 )A<?\x 2 , y 2 ; -A)) (12) 



(where the tilde denotes transposition on the colour indices alone) and we can use 
for and A^ 2 ) the quadridimensional path integral representation of Feynman- 
Schwinger obtaining a corresponding representation for i/ 4 . 

Working on such expressions eventually we arrive at a "second order" nonhomo- 
geneous Bethe-Salpeter equation of the form 

H 4 (x 1: x 2 ; yi, y 2 ) = H 2 [1 \x 1 - yi ) H^\x 2 - y 2 ) - i J d 4 ^ 4 ^ 4 ^ 4 ^ 
H^( Xl - 6) H?\x 2 - 6) 6; m, V2) H A ( m , V2 ; Vl , y 2 ). (13) 

Here H 2 stands for a kind of colour independent one particle propagator and / denotes 
the appropriate kernel obtained as an expansion in the strong coupling constant 
0£ s — f - an d in the string tension a (better in aa 2 , a being a characteristic length, 
typically the radius of the particular bound state). 

We stress that the method by which (13) has been obtained is essentially nonper- 
turbative. 

At the lowest order in a s and era 2 , after factorizing the conservation 5, we can 
write in the momentum space 

I(P1, P2] Pi, P 2 ) = P2] Pi, P 2 ) + iconfCPl, P2] P[, P 2 ) (14) 

(Pi +P 2 =Pi +P2), with 

7 pert = 167r^a s {D pa (Q)qU^ - \<{^Q V - KQ»)<&D P „{Q) + 

+l<($;Qv - KQM d pAQ) + 

+^r^r 2 (S^Q ul - S^MiS;^ - 6t 2 Q, 2 )D pa (Q)} + ... 



and 
with 



(15) 

/conf = jd 3 re tClr J(r,q 1 ,q 2 ), (16) 



J(r, qi, q 2 ) = ^ \Q20 \l 9io ~ Qt + Qw \l Q20 ~ Qt + 
Q10 + Q20 L 



, liolio/ ■ |Qt| . • |Qt| 
+- -(arcsm h arcsm ■ 



|Qt| Q10 Q20 

T\ v q w qiur k . (r 2 u q w q 2u r h 
I - — , 



j + 

r- ak \ q ^ rk - 2(T <£^i + . . . ( i7) 



In (15)— (17) it has been set 



Pl+Pl P2+P 2 n , , ns x 

Ql = o ' 32 = ~ , Q = Pi ~ Pi = P2 - P 2 , (18) 



D pa (Q) denotes the usual gluon free propagator and the center of mass frame is 
understood (qi = —q 2 = q , q^ = (S hh — r h r k )q k ). 
Eqs. (13)-(|PT|) are the main results of our paper. 

Notice, obviously, that, instead of ( |13"D we could have written the homogeneous 
equation 

M*0 = -* / -0jzH 2 ( m P + k)H 2 (r, 2 P - k)I(k, k'; P)^p(k') , (19) 

which is more appropriate for the bound state problem. In this equation rjj = m n ^ m , 
P denotes the total momentum pi+p 2 , k stands for the relative momentum r] 2 pi — r]ip 2 
(qj = r]jP + and in the CM frame q = ), § P {k) is a kind of Bethe-Salpeter 
wave function in the momentum space. 

From (|T9|) by replacing H 2 (jp) with the free propagator 2 ~^ 2 and performing 

an appropriate instantaneous approximation on / [consisting in setting: Qq = 
and q j0 = ^jl^L or PjQ = p '. Q = or ^ _ _ ^^Ll^I _ ^i^±^ an d p Q = 

t;(w' 1 +Wi + w 2 + W2), with Wj = \Jrn 2 j + k 2 , w'j = \Jm 2 + k' 2 ] we can further obtain an 
effective square mass operator for the mesons (in the CM frame P = 0, P = (tub, 0)) 

M 2 = Ml + U (20) 



with M = J m\ + k 2 + Jm?, + k 2 and 



k u k' = j7r^\H /inst k ' k V o / / • 21 

(2^)^ y 2wiw 2 V 2w[w 2 



The quadratic form of Eq. (|2l]|) , obviously derives from the second order character of 
the used formalism. It should be mentioned that for light mesons this form seems 
phenomenologically favoured with respect to the linear one. 
In more usual terms we can also write 

M = M + V (22) 

with 

(k\V |k'> = , * . 1 J inst (k, k') + . . . (23) 



(2vr) 3 4 



WiW 2 W\W 2 



where the dots stand for higher order terms in a s and aa 2 and kinematical factors 

2 

equal to 1 on the energy shell are neglected. In the limit of small ^ the potential V 
as given by (p3|) reproduces the semirelativistic potential of ref. IT On the contrary, 
if one neglects in V the spin dependent terms, one reobtains the hamiltonian of the 
relativistic flux tube model I for quarks without spin with an appropriate ordering 
prescription i. 

Finally we want to mention that a result strictly related to our one, but directly 
in hamiltonian form ([22]), has been obtained by Simonov and collaborators @. 



The following part of the paper is organized in this way. In Sect. 2 we illustrate 
our nonperturbative method on the model case of two spinless particles interacting 
through a scalar field, to which also the usual perturbative derivation applies. In 
Sect. 3 we obtain the mentioned path integral representation of H4 and in Sect. 4 we 
sketch a derivation of Eqs.(13)-(17). Finally in Sect. 5 we discuss how the qq potential 
is reobtained from Eq. (23). 



2. Bethe— Salpeter equation for scalar particles 

Let us consider two scalar "material" fields and 02 interacting through a third 
scalar field A with the coupling ^{(pfA — (p 2 A). Then, after integration over (pi and 
4>2, the full one-particle propagator can be written as 

s , , . \ / / \ 1 \ , A , ^ fVAe iS °^M(A)iA(x,y:A) , x 
G 2 {x-y) = (O|T0(x)0(y)|O) = (iA(x,y;A)) ee J - f VAe is (A )M{A ^ > ( 24 ) 

where A(x,y,A) is the propagator for the particle in the external field A, So{A) is 
the free action for the field A and the determinantal factor M(A) comes from the 
integration on the fields <fi 



= r detg^ + m* T gA) 



1 - 



j=l,2~ det(<9^<9 M + m 2 
~\ £ { T 9 J d 4 xA(x) A F (0) - X -g 2 J d 4 xd A yA(x)A F (x - y)A(y)A F (y - x) + . . .} 

(25) 

(where the upper minus sign refers to the quark , j = 1 and the lower plus sign refers 
to the antiquark (j = 2)) A F denoting the usual scalar particle free propagator. 
The covariant Feynman-Schwinger representation for A(x, y; A) reads 

1 roo rx rs | | j 

A{x, y; A) = -- ^ ds J VzVpexp {1 ^ dr[-p^ + -p^ - -m 2 + - ± gA{z))} 

2 roo rx rs J 

= -- / ds Vzexp{-t dT-[(z 2 + m 2 )TgA(z)]}, (26) 
2 Jo J v Jo 2 



where again the upper sign refers to the quark and the lower to the antiquark; the 
path integrals are understood to be extended over all paths z 11 = z^{t) connecting y 
with x expressed in terms of an arbitrary parameter r with < r < s. In Eq. (|26|) z 
stands for ^f^- and the "functional measures" are assumed to be defined as 

(IT 

vz = (-^) 2 ^V..rfV-i, vp={^f N d'p l ...d i p N . l d i p N 

2-kie 2n 

VzVp = {^-) 4N d A Pl d 4 Zl • • • cPpn-kPzn-kPpn. (27) 
2n 



(e being the lattice time spacing and e — > being understood). 
Replacing Eq.(^) in ( f24"|) we obtain 

G 2 (x-y) = \ Hds fvzexpi— f dr(z 2 + m 2 )}(exp ^ f ' drA(z)), (28) 

2 JO Jy Z JO 2 Jo 

and, if we take simply M(A) = 1 (quenched approximation), 

(exp ^ f drAiz)) = exp ^ f dr f dr'ZU,? - z'), (29) 
2 jo 4 jo Jo 

D F {x) being now the free propagator for field A. 
Similarly for the two-particle propagator we have 

G 4 {xi,x 2 ;y 1 ,y 2 ) = (O|T0 1 (xi)0 2 (^2)0i(?/i)02(?/2)|O) = (G i 2 1 \x 1 ,y 1 ;A)G ( 2 2 \x 2 ,y 2 ;A)) = 
{-f ds 1 ds 2 Vzxl Vz 2 xexp—{ dr 1 (z 2 + ml) + dr 2 (z 2 2 + m 2 2 )} 

2 Jo Jo Jyi Jj/2 2 Jo JO 

(exp ~{g r dr 1 A{z l ) - g f dr 2 A{z 2 )}). (30) 
2 Jo Jo 

and (always for M(A) = 1) 

(exp -{g / dnA{zi)-g I dr 2 A{z 2 )}) = exp ^ — / drj / dTj[D F (zj - z'j)}. 

Zi J JO j ^ 2 * ^ " 

(31) 

In conclusion we have 

^ pOO rCO pX\ PX2 

G 4 (xi,x 2 ;?/i,2/2) = (t:) 2 / ds! / ds 2 / P^i / Vz 2 

2 JO JO Jyi Jy 2 

exp ^[ / 1 dri(i? + ml) - ^ [ 1 dn [ 1 dr[D F ( Zl - z[)] x 
2 Jo 2 jo Jo 



exp — [ / dr 2 (z\ + m\)} - / dr 2 / dr' 2 D F (z 2 - z' 2 )\ x 
2 Jo 2 Jo jo 

exp^i^ 2 -/ dri / dr 2 D F (z 1 — z 2 ). (32) 
4 Jo jo 



stands for Zj(rj), z'j stands for ^(rj)). 
From (32), using the identity 



—ig 2 f Sl f S2 
exp — — / dr ± / dr 2 D F (zi - z 2 ) == 1 - 
4 Jo jo 

- % ^- I dn dr 2 D F (z 1 - z 2 ) exp[— j?- / dr{ / £Z^D P (4 - 4)L (33) 
4 Jo Jo 4 Jo jo 



we can write after some manipulations 

G 4 (x 1 ,x 2 ;yi,y 2 ) = G 2 (x 1 - y 1 )G 2 (x 2 - y 2 ) r(o) / dr i / ds i / rfr 2 / rfs 2 

4 z jo Jti jo Jt 2 

/[' f Xl f X2 f Zl f Z2 
d A Zl d A z 2 Vz x \ Vz 2 V Zl Vz 2 D F ( Zl -z 2 ) 
J J 21 J zo. J vi Jv>. 



'21 -i Z2 Jyi Jy2 

-.2 



exp l -{ dT[[-z'? -m\ + 9 - r dT?D F (z[ - z'()] + P dr' 2 \-z' 2 -ml + 9 - f ° 

Z J T\ Z J T\ J T2 Z Jt'2 

D F (z' 2 - 4')]} x exp ^{ f 1 dr{[-i; 2 - m? + £ f* <D F (4 - 4')] + 
z Jo z Jo 

1 1 2 

+ d7 *H? - ™2 + +^ P d^Dp(4 - 4')] - v dr[ r dT>D F (z[ - 4)} x 
•/() " 2 jo " 1 jo jo 

exp*{^ r dr[ I" dT?D F (z[ - z'{) + £ P P ^'D F (4 - 4) + 

4 Jti JO 4 JY 2 JO 

"T /J ^^(4 -4)}- (34) 

Then, if we replace the last exponential with 1, we obtain immediately the Bethe- 
Salpeter equation 

G 4 (x 1 ,x 2 ,y 1 ,y 2 ) = G 2 (x 1 - yi)G 2 (x 2 - y 2 ) - i J d 4 z 1 J d 4 z 2 J d 4 z[ J d A z 2 

G 2 {x! - Zi)G 2 (x 2 - z 2 )I(z 1 , z 2 ; 4, 4) G 4(4> 4) Vi, V-i) ( 35 ) 

with the ladder approximation kernel 

I(z 1 , z 2 , 4, 4) = 9 2 D F {z 1 - z 2 )5 4 (z 1 - z[)5 4 (z 2 - z' 2 ). (36) 

On the contrary, if we had considered the entire expansion of the last exponential in 
(33), we would have obtained 

I(zi, z 2 , z[,z 2 ) = g 2 D F (z[ - z 2 )S 4 (z 1 - z[)5 4 (z 2 - z' 2 ) - 
-ig A J d%D F ( Zl - z[)G 2 ( Zl - 6)G 2 (£i - 4) 

Df(Zi ~ z 2 )5\z 2 - 4) - ig 4 J d&\ Zl - z[)D F ( Zl - &)G 2 (z 2 - £ 2 )G 2 (£ 2 - 4) 
D F (z 2 - 4) + ig 4 D F (z x - z' 2 )G 2 (zi - z[)D F (z 2 - z[)G 2 (z 2 - z 2 ) + . . . (37) 

(see Fig.l). Finally one can go beyond the quenched approximation and take into 
account additional terms in Eq.(25). This would amount to insert <p\4>i and 2 <^ 2 
loops in all possible ways inside the graph. 

Notice that the final form of the kernel we have obtained is expressed as an 
expansion in the coupling constant g as in the ordinary derivation. However in the 
method described, once we have written eq. (31) and set the last exponential in (34) 
equal to 1, Eqs.(35)-(36) follow exactly. So perturbative expansion appear only at the 



level of successive corrections and not in the basic approximation. This is the reason 
why the method applies even to QCD when we replace (31) with (2). In fact, as we 
have already mentioned, we shall see that in such a case the kernel I is obtained as 
an expansion both in a s and aa 2 . Obviously by Fourier transform of Eq.(35) we may 
pass from this to the more usual momentum counterpart. 

For subsequent developments it is important to mention that we can also obtain 
the B-S equation directly in the momentum space starting from the first line of (26) 
and working in the phase space rather than in the configurational one. The only 
difference in such a case is that we need to make explicit reference to the discrete 
form of the path integral (cf.(27)) and use the discrete counterpart of (33); we refer 
to for details. In the basic approximation we find 

I(p 1 ,p 2 ,p' v p / 2 )=Ag 2 j d^z 2 J / rf 4 ^rf 4 4e-^- fcl )-e-^- fe )- 

D( Zl + z ' 1 - ^1±A) e i(pi-fci)4 e i(pi-fc2)4 (38) 

which would literally correspond to the mid-point discretization prescription (even if 
immaterial in this case). In (37) I{pi,P2',Pi,P 2 ) is t ne ordinary Fourier transform of 
I(z±, z 2 ; z[, z 2 ); introducing the total momentum P = p\ +p 2 and the relative momen- 
tum k = — ^ — pi "^r — £>2 = 7?i Wi — r?2P2, we can also factorize the 5 conservation 

mi+m2 ini+rn,2 I ixrx i^f^i 

function 

I(pi,P2;p' 1 ,p' 2 ) = (27r) 4 5\P-P')i{k,k';P) (39) 

and write 

I{k,k';P) = g 2 D F (k- k'). (40) 



3. The quark— antiquark propagator 

Let us come back to the QCD case. 

The Feynman-Schwinger representation can now be written 

i r 00 is 1 i f 00 

A CT (x, y; A) = -- J dsP xy T xy exp —(—D^D^ - m 2 + -ga^F^) = -- J ds 

[* VzP xy T xy expi I* dr{-)-{m 2 + z 2 ) + gA p {z)z p + 9 -y»F^{z)}, (41) 

Jy JO Z 4 

where P xy and T xy prescribe the ordering along the path of the colour and of the spin 
matrices respectively. 

Furthermore, notice that, as a consequence of a variation in the path z M (r) — > 
z^(t) + 5z^(t) respecting the extreme points, one has 

5 {P xy expig f drz^i^A^z)} = 

J 

= ig jT SS^(z(r))P xy { - F, u (z(r)) exp ig jT dr' z^r')A^z(r'))} (42) 



with 5S^ u {z) = ^{dz^5z u — dz v bz^\ So one can write 

T xy exp(~ J^dra^jJ^-^P^expzg jf' dr> ' z»{t>)A»{z{t') 

= T xy P xy expig J* dT[z»(T)A^z(T)) + ^ F^{z{r))] (43) 
and Eq. fl4l|) can be rewritten as 
A ff (x,y;A) = -~ f S dr [ X VzP xy T xy S^exp i f dr[-\{m 2 + z 2 ) + ig^ A^z)} (44) 

Z JO Jy JO £ 

with 

So=exp[-i [dra^— 5 —-\ (45) 



A Jo 5S^(z) 

In (44) it is understood that z M (r) has to be put equal to z^(t) after the action of 
Sq. Alternatively, it is convenient to write z = z + (, assume that S$ acts on C(r) 
with SS^(z) = \{dz^5C - dz u 5(v) and set eventually C = 0. 
Replacing (44) in (12) we obtain 

[ poo poo rxi rX2 

H 4 (xi,x 2 ;yi,y2) = (-) 2 / dsi / ds 2 / Vz 1 j Vz 2 T Xiyi T X2 y 2 SQ 1 S S2 

l JO JO Jyi Jy 2 

exp (^){ f 1 dn(m? + i') + f * dr 2 (m^ + i 2 2 )} 
I Jo Jo 

i(TrP r exp(^){ jf (46) 

where now z = ^ = z 3 - + on 1^ and T 2 , z = z on the end lines X\X 2 and ?/ 2 2/i an d 
again the final limit Q — > is understood. 

Then, let us try to be more explicit concerning Eqs. (2) and (3). For the first 
term in (2) we have, at the lowest order of perturbation theory, 



i(ln W) pert = zln(~TrPexpz# f dz^A^z))^ = -g / dr x \ dr 2 D ilv {z 1 - z 2 )z1 

3 Jr 6 Jo Jo 

-?9 2 / dn / diiD^zi - z'Jzfz'" - -g 2 / rfr 2 / dT' 2 D^(z 2 - z' 2 )z^z 2 + . . 
i Jo Jo i Jo Jo 



On the other side, for the second term in general we have to write 

x m _ u^(t, A) being the equation of the minimal surface with contour V. Let us assume 
that for fixed t and for A varying from to 1, « M (A, t) describes a line connecting a 
point on the quark world line 1^ with one on the antiquark world line T 2 , 

u"(M) = zZ(n(t)), u"(0,t) = z$(r 2 (t)). (49) 



Obviously (^) is invariant under reparametrization, so a priori the parameter t could 
be everything. In particular, however, if I\ and T 2 never go backwards in time, t 
can be choosen as the ordinary time, as in (3), u°(s,t) = t. Then Ti(£) and T2(t) are 
specified by the equation 

4(n) = 4(r 2 ) = t. (50) 

We shall set 

Obviously L cannot depend only on on the extremal points Z\{t\) and z 2 (t 2 ), but has 
to depend even on the shape of the world lines, at least in a neighbourhood of such 
points. So, we can think of it as a function of z±, z 2 and of all their derivatives in t\ 
and r 2 , L = L(z±, z 2 , Zi, i 2 , . . .). Finally we can write (48) as (recall Zj 



2° / d^^-^L^i^.iijia,...) = J dn J dT 2 5(zi-z%)ziZ%L(z u z 2 ,z 1 ,Z2, 

(52) 

In principle this expression can be considered a good approximation even if the world 
lines contain pieces going backwards in time. In fact, in such a case if we fix e.g. T\, 
( |5H| ) has more than one solution in r 2 and, if Y\ and T 2 are not too much irregular 
in space (otherwise S m i n is large and the weight of the loop is small), the minimal 
surface can be reconstructed as the algebraic sum of various pieces of surface. 
In the straight line approximation we must choose 

u\\t)=t ■ u\\t) = \z k (r 1 (t)) + (l-X)z k (r 2 (t)) (53) 

and we have 

rl 1 

zJigL = ff|zi - z 2 | / d\{zl zl - (Az 1T i; 2 o + (1 - A)z 2T i 10 ) 2 }2 (54) 
Jo 

which, introduced in (52) becomes equivalent to Eq.(3). The important point is 
that (52) with (54) has the same general form as (47). However we stress that the 
approximation (53) must be performed only after that the application of the operators 
iSq 1 and Sq 2 has been performed. 

Substituting fl4"7p and (52) in (HBp we obtain at the lowest order 



[ roo rco rxi rX2 

H i {x u x 2l y u y 2 ) = (-) 2 / ds 1 ds 2 Vz 1 Vz 2 r T Xiyi r T X2y2 S s a 1 S { 

L JO JO Jyi Jy2 

1 r Si 1 

eX P^{ - o / dl ~l( m l + z l) - o / dT 2( m l + Z 2 ) + 

2 Jo 2 Jo 



-I s -fl --VI 
Z 2 ) Z 2 Z 2 



+ o9 2 / dn / dr^D^Zx - + -g 2 / dr 2 / dT' 2 D^ u (z 2 

6 Jo Jo 3 Jo Jo 

- dr x \ dr 2 E(z u z 2: zi,z 2 ,...)}, (55) 
Jo Jo 



where we have set 

E(zi, z 2 , 21, z 2 . . .) = ^g 2 D flu (z 1 - z 2 )iii£ + 

+cr<5(>io - Z2o)z 10 Z2oL(z u z 2 , zi, z 2 . . .)■ (56) 

Now, let us denote the quantity in curly bracket in (55) by $ and perform a 
Legendre transformation by introducing the momenta p J(U = — 4™ ( where the various 

quantities assumed to be treated as independent) 



2 Q-M 



Pui = z^ + ^g J d^D^zx-z'xYzx + J dr. 

= ^2 + -£ y o dr 2 D^(z 2 - z 2 )z 2 + dr x ^ . (57) 

Eq.(|57j) cannot be inverted in a closed form with respect to z\ and z 2 , however, we 
can do this by an expansion in a s = ^- and aa 2 and at the lowest order we have 

•/* * 4 2f Sl j'n / f 2 A t dEjzuz'^pup^...) 

Z\ = y o dr^D^Z! -z x )v x - dr 2 ^ + ... 



4 2 >n I- -'\ n, f S1 A ,dE(z 1 ,Z 2 , P , 1 ,p 2 ... 



z 2 = P2~^9 J Q dT 2 D^ v {z2- z' 2 )p 2 - dr[ ' ^ ' — + ...(58) 

with 



l>", = Pf + Cf (59) 

In conclusion we find (up to a determinantal factor that in this approximation can 
be set equal to 1) 

[ poo poo pX\ pX2 

H 4 (x 1 ,x 2 ,yi,y 2 ) = (~) 2 / ds x ds 2 Vz 1 Vp 1 Vz 2 Vp 2 T Xiyi T X2y2 

Z JO JO J yi Jy 2 

Usi rS2 rsi rS2 ■> 

dr 1 K 1 + J^ dr 2 K 2 - J dr x dr 2 E{z x , z 2 , pi, p 2 , ...) + ... j 

(60) 

where 

1 ^ pSj 

K 3 = ~Pj ■ + 2^~ + 3 g2 Jo 3 dT j D ^(zj ~ tyPjPj' + • • • ( 61 ) 

includes the self-interacting term. Notice that now in Sq j it must be understood 
SS^) = Idr^dq - p»5$) + .... 

4. The Bethe— Salpeter equation in QCD 

From Eq. (|60|) we proceed along the same line followed in Sect. 2. with reference 
to Eq.(32). 



Using 



exp / dn / dT 2 E(z 1 ,z 2 ,p u p 2 . . .) = 
Jo Jo 

= 1+ dr 1 E(z u z 2 ,p 1 ,p 2 ...)exp dr E(z[,z 2 ,p l ,p 2 ,...) 
Jo Jo 

corresponding to (33), we have 



(62) 



H 4 (x 1 ,x 2 ; 2/1,2/2) = (^) 2 



T Xl y 1 T X2 y 2 S 1 S 2 ^ expz 



ds 



dr\K^ 



dr 2 K 2 \ 

S2 



J-2 



VzxDpx / Vz 2 Vp 2 



driKi + / dr 2 K 2 — / g?t( / dr 2 E(z[, z 2 ,p'i 
Jo Jo Jo 



dn / dr 2 E(zi,z 2 ,p h p 2 . . .) 



(63) 



L ' ' ' " 1P11 ■ ■ ■) 



To obtain from this an equation analogous to (34) we need to commute S^Sq 2 with 
E. To this aim using the method of Ref. u and bearing in mind (|56|) , (|5l| ) and (|55|) , 
we find first 



5S^(z 1 ) Jo 



rsi rs 2 

/ dr[ / dT l 2 E{z' l ,z' 2 ,p' ll p 2 , . . .) 

JO JO 



+a5(z w - z 20 ) 



(pio - pi)(zi - Z2) 2 + (pi ■ (zi - z' 2 )y 



+ 



(64) 



and a similar result, with a minus sign in front, for the derivative SSfl t( z / ) ■ Furthermore 



SS^iz^SP^z'^ Jo 
but 

S 2 

8S^{zi)8S^{z 2 ) 



dr'l 



S2 11 

dr i E= 5S^(z 2 )5Sp°(z' 2 ) 



- P dr'{ r dr'lE = 0, (65) 

9 JO JO 



(66) 



Then, taking into account the relation 

00 1 



e A Be -A = ^ [A [A [A,B] ...]], 
^0 n! 



(67) 



we have 



T 2 +e\-l 

T2—E J 



ds' 2 a^ 2 



7t) 



J dn J dT 2 E(z 1 ,z 2 ,pi,p 2 ,...) = R{zi,z 2 ,p h p 2 ) 



(6* 



with 

R = -Rpert + -Rconf (69) 

4 f 

~\<^d lu - SZd^Drfa - ^)p 2 - \<(5;d 2v - 5ld 2 ,)D pa ( Zl - z 2 )p{ 
+ * rvr2(( ^ _ d P idlpi )(5^d 2u2 - d° V2 d 2li2 )D pa ( Zl - z 2 )} (70) 



and 



-Rconf = Cr^io - 2 20 )||Z1 - Z 2 | ds^pl Q pl - [sp 1T p 20 + (1 - s)p 2T Plo] 2 

1 fiiyPlu(zi p — Z 2p ) — Pl p (zi u — Z 2v ) 1 nvP2v{Z\ p — Z 2p ) — p 2p (zi u — Z 2v ) 

— -,V2qO\ 1 V TP1O02 



|zi -Z2IVP10 -Pit 4 " l z i - z 2|vt4-pIt 



(71) 



Notice that in (68) we have eventually supressed reference to the higher order deriva- 
tives in z\, z 2 and this corresponds to the adoption of the straight line approximation. 
Finally setting 

H 2 (x-y) = — / ds VzVpT xy S^expi drK (72) 

2 JO Jy JO 

we can write Eq.(18) as 

H 4 (x 1 ,x 2 ; yi, V2) = H 2 (xi - y{)H 2 (x 2 - y 2 ) + 

2 roc rco rx\ rX2 rs\ rS2 

— \ ds x \ ds 2 / Vz{Dp x \ Vz 2 Vp 2 T Xiyi T X2y2 / dr x \ dr 2 R{zi, z 2 ,p 1 ,p 2 ) 

4 JO JO Jyi Jy 2 JO JO 

S^S^ exp i{ £ dr[K[ + jH dr' 2 K' 2 - 1 jT dr[ jH dT>E(z[, z' 2lPl ,p 2 ...)}• (73) 

At this point, as mentioned at the end of Sect. 2, to go ahead it is necessary to 
take explicitly into account the discrete form of (|73|). If, in particular, we set 

2 ' 
(74) 

(as required by a gauge invariant definition of the integral on the gluon field) by 
appropriate manipulations we eventually obtain (13) with (cf.(38)) 

(27v) i 5(p 1 +p 2 -p' 1 -p' 2 )i( Pl ,p 2] p' 1 ,p' 2 ) = -U J d%d% J d%d% 
dki d 4 k 2 jfo-fcjfr-Hfa-fcajfa p/ 6 + rjx C2+V2 



Pexp [ig I dz»A p {z)} =P'[[U{z n ,z n - 1 ) = Pex V igJ2K ~ 4-i)M 



(75) 



Then, using (69)-(71) and performing explicitely the integrals we find Eqs. (14)-(18). 

Obviously the homogeneous Eq. (19) follow from (13) using the usual decompo- 
sition of the quark-antiquark propagator in terms of Bethe-Salpeter wave functions 

GJk, k', P) = V + regu i ar terms. (76) 

b P 2 -m B 



5. The semirelativistic potential 



From (23), expanding in ^ and passing to the coordinate representation, we obtain 
the following potential 

4 a, 4 - 1 - h - k 



V 



a. 



h ar H — 

3 r 3 mi?7T-2 2r 



{ 1 {5 hk + ^r )qhqk] 



w 



4 



1 ol\ • r 1 a 2 • i\ 

2m 2 r 3 



1 



3 " x 2mi r 3 



3 2mi7«2 



>i + 02) • (r x q) + 



3 mim 2 

a , 1 1 
- — + — 

6 mf 



[ 3(ai ' r) ^ 2 - r) - ^] + ^^^1 • a 2 )5 3 (r) 



3 mim 2 3 



m 2 mim 2 
m 2 ' / 



-){QT r }w 



+ =3J ■ 17 x qj - -1 



ai 1 ai • i" 1 a 2 • r 



(77) 



2 mi r m 2 r 

where now q stands for the momentum operator and the symbol { }w stands for 
the Weyl ordering prescription. Notice the nonhermitian terms in the Dirac matrices 
«i and a 2 . Such terms can be eliminated by performing the Foldy-Wouthuysen 
transformation with the nonhermitian generator 



-a x ■ q - 



2mi 2m 2 
and we end up with the semirelativistic potential 
4a 



a 2 • q 



3 r ^ 2/;; 



+ 
+ 
+ 



^\^ {6 ^ + f h f k )q h q k \ 
\rri2 13 r J 



w 



Et- 2 { (7r 4}w + 

• =1 U»lj 



1 



6mim 2 



{arq^jw + 



2 

mf 



ml 



4 a, 



3 r 

^Si • (r x q) + -^S 2 • (r x q) 



4 a, _ 
3r 3 
1 4 01 

-%.(rxq) + S r (rxq)] + 



m 



rrin 



+ 



mim 2 3 r 3 
1 4 



(78) 



mim 2 3 



■-a, < — 



1 



-(Si-r)(S 2 -r)-Si-S 2 



+ y5 3 (r)S! ■ S 2 } (79) 



which coincides with the Wilson loop one I. 
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